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Abstract 

We study correlation functions of Wilson loops and local operators in a sub- 
sector of N = 4 SYM which preserves two supercharges. Localization arguments 
allow to map the problem to a calculation in bosonic two-dimensional Yang-Mills 
theory. In turn, this can be reduced to computing correlators in certain Gaussian 
multi-matrix models. We focus on the correlation function of a Wilson loop and 
two local operators, and solve the corresponding three-matrix model exactly in 
the planar limit. We compare the strong coupling behavior to string theory in 
AdS$ x 5 5 , finding precise agreement. We pay particular attention to the case 
in which the local operators have large R-charge J ~ vA at strong coupling. 
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1 Introduction 



Exact results in non-abelian gauge theories are rare and clearly of great importance. In 
supersymmetric gauge theories, the powerful technique of localization allows sometimes 
for such exact results for certain observables preserving some fermionic symmetries of 
the theory. For example the exact expressions conjectured in [l|, 0] for the circular 
1/2-BPS Wilson loop in J\f = 4 SYM were proved using localization in jsj], as well as 
extendend to a large class of M = 2 theories. 

In M = 4 SYM, a wide generalization of the 1/2-BPS circle to lower supersymmetric 
Wilson loops of arbitrary shapes was introduced in jjj^] and then classified in Q. 
An interesting subfamily of that construction consists of operators supported on any 
loop on a two-sphere S 2 embedded into the M 4 spacetime. Generically, these Wilson 
loops are 1/8-BPS, and it was conjectured in [4j-l6( that their quantum correlators are 
exactly captured by a purely perturbative calculation in bosonic 2d Yang-Mills. The 
conjecture was later strongly supported by the localization calculation in [8j, where 
it was shown that the path-integral with insertions of those loop operators localizes 



on a 2d gauge theory closely related to the Hitchin/Higgs-Yang Mills system 
which can be seen to be perturbatively equivalent to ordinary bosonic 2d Yang-Millso 
The calculation in the 2d theory can be then mapped to certain Gaussian multi-matrix 
models, which allow for an exact evaluation of the correlators. Several checks of the 
relation to 2d YM have appeared 001 In particular, recently the localization 



result for the expectation value of a Wilson loop in this family was used in [18 



sec 



also [19|) to derive an exact expression for the low-angle limit of the cusp anomalous 



dimension. This has been checked u sing integrability up to three loops in 20] (see also 



21]) and analytically to all loops in 22], providing a first link between localization and 
integrability results. 

The calculation in Q also suggested that localization applies in fact not only to 
the Wilson loops, but to a larger sector of operators that are annihilated by the same 



as 



supercharge. This include certain chiral primary operators inserted on the S 2 [23 
well as 't Hooft loops linked with the S 2 24|. The correlation function of a Wilson loop 
and a local operator in this sector was computed in [23j], giving support and generalizing 
the original conjecture of 25] for the correlator of a 1/2-BPS Wilson loop and a chiral 



orimary (see also 26] for the study of the large R-charge limit of this correlator). In 



24 



27] . the exact results implied by the relation to 2d were also used to obtain some 



2 In [8j the one-loop determinant for fluctuations around the localization locus was not computed. 
The conjecture of [4|-|6j follows if one assumes that the determinant is just unity. Hence the conjecture 
for the 1/8-BPS loops is not yet proved at the same level of rigour as 

3 A11 calculations in the literature agree with the original conjecture except for a certain disagree- 
ment found in 16]. It would be good to clarify the nature of that disagreement. 
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explicit tests of the S-duality symmetry of the M = 4 SYM theory. 

In this paper, we continue the study of this supersymmetric subsector and concen- 
trate on mixed correlation functions of Wilson loops and local operators. In particular, 
we focus on the limit in which the local operators have large R-charges J ~ y/\ in 
the strong coupling regime A ^> 1. A motivation to look at this problem is the re- 
cent progress in computing correlation functions of operators with large charges using 
semiclassical strings in AdS (see e.g. 28l438j and references therein). We hope that 
studying similar correlators in a subsector where exact results are possible may pro- 
vide detailed tests of the gauge/string dictionary and serve as a useful benchmark for 
various computations involving more general non-protected operators^ 

As mentioned above, the localization to 2d allows to map the calculation to certain 
multi-matrix models. Atthree-point level, we have to solve a Gaussian three-matrix 
model. Since we could not find explicit results available in the literaturejf] here we 
present in detail a derivation of the exact planar resolvent for the most general Gaussian 
three-matrix model. The result can be applied to all possible three-point correlators of 
Wilson loops and local operators, but in this paper we focus our attention to the case of 
(W O j x O j 2 ) . After deriving an exact prediction for this correlator, we study its strong 
coupling behavior in the regimes Ji, J 2 <C 1> Ji ~ \/A ^> J 2 and Jj, J 2 ~ vA- In the 
first case, in particular, we observe a factorization at leading order at strong coupling 
which is completely analogous to the one seen in [44| for four-point functions of two 
"heavy" and two "light" operators. For J 1( J 2 ~ y/\ we derive a prediction for the 
exponential behavior of the correlator which should be matched against a semiclassical 
string solution with two spikes which end on two boundary points. 



On the string theory side, we first derive a generalization of the solution of 26 



describing (WOj) for J ~ \/\ ^> 1 to the case of the 1/4-BPS circular loop. The 
corresponding string lies in Ad S3 x S 3 and preserves less supersymmetries than the 
solution of [26]. Computing its area we find precise agreement with the two-matrix 
model derived from the 2d YM description. We then use this solution to obtain the 
string prediction for (WO^O^) at J\ ~ vA 3> J 2 , again observing perfect agreeement 
with the exact solution of the three-matrix model in the appropriate limit. In this 
paper we do not find the string solution which should describe (WO^Oj^) in the 



4 Some examples of (non-supersymmetric) 

correlation functions of Wilson loop and local operators in the limit of large charges were recently 
studied also in 39| 40|. 

5 There is a vast literature on multi-matrix models, see e.g. 41 1 42 1 43[. The class of models that 
can in general be solved exactly are those of open chain type, where the interactions among the n 
matrices are of the form XiXi+i, i = 1, n. In our case we also encounter matrix models with closed 
chain interactions which are not solved in general. However, since our models are Gaussian, a solution 
is still possible, as we show. 
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regime Ji, J2 ~ VX We derive however a set of first order differential equations which 
follow from supersymmetry of the system and which should be satisfied by all string 
solutions in AdS 3 x S 3 dual to the general (WW • ■ ■ 00 ■ ■ ■) correlators in our sector. 
Hopefully, these equations will be useful to find or characterize new explicit solutions. 

The paper is organized as follows. In Section 2, we review the Wilson loops and 
local operators of interest as well as the dictionary relating them to the 2d theory. 
In Section 3 we show how to derive from 2d YM the Gaussian multi-matrix models 
capturing mixed correlators of Wilson and local operators (in particular (WOO)). In 
Section 4 we solve the general Gaussian three- matrix model in the planar limit, obtain 
an exact prediction for (WOO) and study its strong coupling limit. In Section 5 we 
derive the string solution dual to (Wi/^bpsOj) at J ~ y/X ^> 1, compare its area to 
the localization prediction and use it to compute (WOO) in the limit of one "heavy" 
and one "light" local operator. Finally, we study the supersymmetry constraints on 
the string solutions in our sector and derive a set of first order differential equations 
that the solutions should obey. 



2 Review of supersymmetric subsector and relation 
to 2d 

In our conventions, the M = 4 SYM action on R 4 with the standard flat metric is 

SsYM = 2 / d ^ X ( o tr 

+ trD^ A D^ A + ...) , (2.1) 
9ym J \ z / 

where /1 = 1, . . . , 4 are space-time indices and A = 1, ... ,6 are SO(6)r indices. The 
covariant derivative is D = d + A, the curvature is F^ v = [D^, D u ], all fields take value 
in the Lie algebra of the gauge group U (N) and represented by anti-Hermitian matrices 
A^ = A^Ta, $ = $ a T a . The anti-Hermitian generators satisfy tr T a T& = — |5 a &, hence 
the action may be also written as 

Ssym = 2jjr- J {^f; u f; u + D^\D^ A + ...). (2.2) 

The Wilson loops we study in this paper are the 1/8-BPS operators constructed in 
0-0]. They are supported on arbitrary closed curves on a S 2 inside M 4 which we may 
define in Cartesian coordinates as 

3 

x 4 = 0, $>, 2 =r 2 . (2.3) 

i=l 
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The radius r of the two-sphere is arbitrary, and we will set it henceforth to 1 for 
simplicity (the radius dependence is easily reintroduced if desired). The 1/8-BPS 
Wilson loops couple to three of the six scalars, $j, i — 1, 2, 3, and for any loop C G S 2 , 
they are given by 



W R (C) = -J- tr R Pexp <f(Ai + ie ijk ^iX k )dx ] , (2.4) 

«R Jc 

where dn denotes the dimension of the representation R. In other words, given a loop 
defined by x(s),x 2 = 1, the operator couples to the combination of scalars (x x x) ■ $. 
For arbitrary curve, these operators preserve four supercharges^] Supersymmetry can 
be enhanced for special shapes. For example, the well-known 1/2-BPS circular Wilson 
loop is obtained by taking C to be an equator of S 2 . Circles of arbitrary radius along 



latitudes of S are 1/4-BPS and they coincide with the 1/4-BPS Wilson loops of [45 
Note that since the four supercharges preserved by the loops do not depend on the 
contour, an arbitrary collection of Wilson loops on S 2 is also 1/8-BPS. 



As shown in [23[], see also [8j, it is possible to add an arbitrary number of local 
operators on the same S 2 while still preserving two supercharges. The relevant local 
operators are the following 

Oj(x) =tr (x i $ i + i$ 4 ) J x l e S 2 ,i = 1,2,3 (2.5) 

Note that these can be viewed as ordinary chiral primaries inserted at a specific point, 
where the orientation in the scalar space is correlated with the position of the operator. 
It is easy to see from the definition that the two-point function of these operators is 
position independent. In the planar limit we have 

(Oj{x)Ojt{x')) = J (j^r^) Sjj>- (2.6) 



In fact, (12.51) are a special case of the superprotected operators introduced in 46], where 
it was shown that all the n-point functions (Oj 1 (xi)Oj 2 (x2) ■ ■ -Oj n (x n )) are position 
independent and moreover tree-level exact. 

The system of any number of these local operators on S 2 preserves four supercharges 



46] 23]. When the Wilson loops (\2A\i are also present, the combined system is invariant 



under two supercharges [23]. While all n-point functions of the Oj's are protected, 
mixed correlation functions of Wilson loops and local operators can have a non-trivial 
coupling dependence and will be the focus of this paper. 



6 In flat space, these supercharges are linear combinations of Poincare and superconformal super- 
symmetries. 
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Figure 1: Supersymmetric Wilson loops and local operators on S 2 . The 
system of any number of loops of arbitrary shape and any number of 
local operators preserves 2 supercharges. Correlation functions of these 
operators are mapped by localization arguments to a multi-matrix model 
computation. 



A linear combination of the two supercharges preserved bv the system is precisely 
the fermionic charge used in the localization calculation of [8], which therefore also 
applies to mixed correlators of Wilson loops and local operators. Hence, extending the 
conjecture of 0-0], it was proposed in [23] that correlators of any number of Wilson 
loops and local operators can be computed exactly by the bosonic two-dimensional 
Yang-Mills theory on S 2 with action 



S- 



YM 2 



2& 



(2.7) 



with the following 4d/2d map 



W R {C) ti R Yj 



o 



tr ( i *2d 



2 
92d 



9m 
2tt 



(2- 



where we used tilde to denote the two-dimensional fields. Under this map, on the 
2d YM side we should extract only the perturbative, or zero-instanton, contribution. 
In fact, non-trivial instantons on the 2d side also have an interpretation in 4d: they 
correspond to turning on 1/2-BPS 't Hooft loop operators along a S 1 linked with the 
S 2 24j (see also [27)). In this paper we will not consider 't Hooft loops and focus on 
Wilson loops and local o per ators. One way to extract perturbative contribution of the 
exact 2d YM quantities 47H5l| is to decompose them into a sum of instanton sectors 
and then pick the term with zero instanton number 52H54j . Another way is to simply 
sum up the 2d perturbative expansion around the trivial vacuum A = in a gauge in 
which the theory becomes free. Either way, the end result is that the 4d correlators 
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can be eventually mapped to certain Gaussian multi-matrix models 

(W RL {C 1 )W Ra (C 2 ) ■ ■■Oj 1 (x 1 )O j2 (x 2 ) ■■■) Ad (2.9) 
= J [dX 1 ] [dX 2 ] ■ ■ ■ [dYj] [dY 2 ] ■ ■ ■ tr Rl e Xl tr R2 e X2 ■ ■ ■ tr Y/ 1 tr Y 2 h ■ ■ ■ e ~ Sm - m - [x;Y] 

where the matrix model action S m _ m _[X, Y] is a quadratic form in X a , Y a whose coef- 
ficients depend on the areas singled out by the Wilson loops and the topology of the 
system (i.e. in which regions the local operators sit, but not on their precise position). 
Note that we did not assume a large N limit here. Localization, and hence the matrix 
model description, should apply at any finite N . 

As a special case, the Wilson loop expectation value is given by the 1-matrix model 
(W R (C)) M = iy [dX] ti R e x e~^t^ trX2 (2.10) 

where A\,A 2 are the areas singled out by the Wilson loop and A = A% + A 2 = Air. 
In particular for the 1/2-BPS circular loop A\ = A 2 = A/2 and the correct Gaussian 
matrix model [l| j^J is reproduced. 

Another example previously studied is the Wilson loop/local operator correlator, 
which using the map to 2d YM can be shown to be given by |2c 



AC - A tr Mly2 2i XY \ 

(W R (C)0j(x)) M = - J[dX][dY] ti R e x tiY J e ^ Ua " 2 >, (2.11) 

where we have assumed that the local operator is in the region of area A\. This result 
in particular reproduces and generalizes the conjecture of 25[for the exact correlator 
of the 1/2-BPS circular loop and a chiral primary (see also 55] for the generalization 
to the 1/4-BPS circle) fl The correlator of two Wilson loops on S was also studied, 



sec 



14 15 



Note that the position independence and tree-level exactness 46( of the correlation 
functions of local operators (12. 5 j) can also be easily seen from the point of view of the 
2d theoryjfl The explicit multi-matrix model which computes the correlator in this case 
is given by 

(0 Jl O j2 ...0 Jn ) = ±J[dX 1 \[dX 2 \...[dX n \ c -^K^G^i*-> a -£^*2) x 



7 Conformal symmetry fixes the correlation functions of a circular Wilson loop and a scalar primary 
operator up to an undetermined function of the coupling, see e.g. 3^, 5(|. Hence the exact (WO) 



correlator on the sphere is enough to determine the correlator for arbitrary position of the local 
operator. 

8 One may derive for example a Ward identity by acting with the differential on the correlation 
functions of tr(i * 2 d F) J operators to show their position independence. Tree level exactness follows 
from the fact that 2d YM becomes Gaussian in an appropriate gauge. 
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x tr tr X? ■ ■ - tr X J n n . (2.12) 

The matrix-model action is chosen to reproduce the tree-level propagator between the 
local operators, and is such that the propagators from an operator to itself are set to 
zero (the operators are understood to be normal-ordered). Sometimes in this paper we 
will also consider operators normalized in the same way as ordinary chiral primaries 
with unit two-point function, i.ej^l 

Oj= (-i^ -Ltr^ + i&) J = AfjO j. (2.13) 

Let us mention that the localization arguments reviewed above should also apply 
when the local operators are inserted along the loop, i.e. for gauge invariant operators 
of the form 

tr{(4^(xi) + i^(x l )) Jl W[x 1 ,x 2 }(x i 2 ^ t (x 2 ) + i^\x 2 )) j2 W[x 2 ,x 1 }} 

rxrr it, [ X2 (Ai+iE iik $>iX k )dxi (2-14) 

W[xi, x 2 \ = Pe Jx i v 3 ' 

and analogous operators with arbitrary number of insertions along the loop. These 
should map to the corresponding operators in 2d YM according to (12. 8p . We leave a 
detailed study of these operators for the future. 



3 Mult i- matrix models from 2d YM 

In this section we show how to derive from 2d YM the Gaussian multi-matrix models 
computing the correlation functions of Wilson loops and local operators. First, we will 
use perturbation theory in the light-cone gauge to obtain the 3-matrix models for the 
explicit example of a triple trace correlator of a 1/8 BPS Wilson loop and two local 
operators on S 2 . Then, we will present a simple formula giving the multi-matrix model 
for the most general multi-point correlators, and show how to derive it from the known 



exact solution of 2d Yang-Mills 47H51 . 



3.1 Light-cone perturbatione theory 



As done in [23(, a simple way to derive the matrix model is to look at the pertur- 
bative Feynman diagram expansion in the A z = gauge, where z, z denote complex 



3 A chiral primary AO tr(u • §) J with u 2 = and u ■ u* = 1 has unit normalized 2-point function in 

the planar limit if Afj = 2 J / 2 ^— ^Tf- The factor of i in the normalization factor is due to our 
convention that the gauge group generators are anti-hermitian. 
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coordinates on S 2 with metric (the radius is set to 1) 

, 9 Adzdz 
ds 2 



'l + zzf 



In the A s = gauge there are no interactions and the 2d YM action becomes simply!^! 

Sym 2 = 4" / <Pz y /gg' s g* s tr(d g A t d g A t ) . (3.1) 

92d J 

We use notations d 2 z = dz A dz = 2idx A dy for 2 = x + zy, and = —igg z , so that 
d 2 z^/g is the conventional volume form on S* 2 normalized as 

y d 2 Zy /g = 4tt. (3.2) 



The gauge field propagator is 



Using this, one gets the propagator for the field strength i * 2 d F = — 1(1 + zz) 2 dzA z 
(i *2 d Fj(z) i * M Ff(w)) = -5}5 k 3 (j| - ^(1 + z^) 2 5 2 (z - «,)) . (3.4) 

For convenience, let us record here also the propagator between the field strength and 
the gauge field 

(i * 2d F*(z) (A z )f(w)) = -fW * _ 1 + ZnJ . (3.5) 

There are two distinct topologies for the correlator of a Wilson loop and two local 
operators on S 2 : one where the operators are on opposite hemispheres compared to 
the loop and the other where the operators are on the same hemisphere, see Figure [21 

Let us start with the case in which the two local operators are on opposite hemi- 
spheres compared to the Wilson loop. Using the area preserving invariance of 2d YM 
and the position independence, we can always choose the loop to be a circle at some 
latitude angle, and place one local operator, say Oj x , at the north pole and the other, 
Oj 2 , at the south pole. Let us parameterize the loop as z(r) = r e lT , where r = tan | 



10 Recall that we use conventions in which the gauge field is anti-hermitian, as in 23]. Also, in this 
section we will omit tilde's on the 2d fields, since confusion with 4d fields will not arise. 
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(a) (b) 

Figure 2: The two distinct topologies for the correlator (WOO). Fig. (a) 
is mapped to an open chain 3-matrix model, while (b) to a closed chain 
one. 



and #o is the latitude angle. It is then easy to see that the propagators in the 
gauge are all constants (i = d T z): 

<*(4.)JW*(4.)f(*)> = -*#7T&* = - 6 > 6 > 9ii 



2tt (1 + rlY ' ' 8ir 2 A 

{i * M F<(0) HA,){ w> = -a^j^j = -A'^f T 
<i » M Fj(oo) i . M F*(0)> = -6<6*§i . 



(3.6) 



Since the propagators are all constants, the sum of Feynman diagrams is obviously 
given by a matrix model. If we assign matrices X%, X 2 , X 3 to respectively Oj 11 W, Oj 2 , 

Iym = 



the no n- vanishing matrix propagators are then, using the relation g\ d = —2/Ag\ 



2 -A 2 1 



;(Xi)}(X 3 )f ) = 3y M ^^'? = Twisty** 



(3.7) 



where we have introduced the shorthand notation A & to denote the propagator from a 
to b matrix, a, b = 1, 2, 3. Inverting the propagator, one gets the matrix model action 

Sowo = 2gJ^ tT (t 2 x * ~ a\a X " + T X * ~ T 2 XlX2 + T X2Xs ) (3 - 8) 
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and the correlator is given by 



{Oj x W r (C)Oj 2 ) = i / [dX x ][dX 2 ][dXz\ trX} 1 tr R e Xa tr X 3 Jz e~ s °* 



(3.9) 



We see that only adjacent matrices interact. This is known in the literature as an open 
chain multi-matrix model. The same structure arises for the correlator of three Wilson 
loops with the topology of three latitudes on S 2 , see [l4j[l5|- In that case, the explicit 
three-matrix model is given by 



-Swww 

j 



(Wr^Wr^Wr^Cs)) = ^J [dXi][dX a ][dX 3 ] tr Rl e Xl tr R , 2 e x Hr R3 e Xs e 
Swww = (4-Xi + ~ X 2 ) 2 + ^-{X 2 - X 3 ) 2 + ±-Xi) . (3. 10) 

Z 9YM \ A l A 12 A 23 A 3 
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where Ax, v4 12 , v4 23 , A 3 are the areas of the regions singled out by the three Wilson 
loops. The generalization to the case of an y n umber of Wilson loops with the topology 
of latitudes on S 2 is straightforward 14J 15[ and gives an open chain multi-matrix 
model (see also next section for the derivation of the most general correlator). 

When the local operators are on the same hemisphere, a similar derivation goes 
through. Assuming that both local operators described by X\ and X 3 are to the north 
of the Wilson loop described by X 2 , the only change compared to (I3.7P is that now 

<(* a )5(*s)?> = wrM^m - 4^ 

i 1 (3.11) 



Then the matrix model action is 



S oow = 2 ,f , v tr f-A 2 X 2 + 5 - A 2 X 3 2 - 2^^ - 2*X 2 X 3 + 2AX 3 Xx 

(3.12) 

In this case all pair of matrices interact, and we may refer to this model as a closed 
chain. 



3.2 General multi-matrix model from 2d YM exact solution 



In 14j . a simple way to derive the multi- matrix model giving the correlator of several 
Wilson loops was proposed. The argument was based on the observation that 



48 



51| 2d YM on a Riemann surface £ localizes on classical configurations solving 
0. Also, two dimensional Yang-Mills theory essentially reduces to the abelian 



49 

d A *F 

theory 48j, |49j, so that one can take *F to be a piecewise constant function on E with 
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jumps supported on the Wilson loops. The value of *F in each region is then related 
via Stokes' theorem to the integral of the gauge field § c A = iXi along each loop, where 
the matrices are the variables in the multi-matrix model. Let us review how this 
works in the simplest case of a single Wilson loop expectation value on S 2 . We divide 
the sphere in two regions S 1; E 2 of areas Ai, A 2 . The field strength is taken to be 

* F = i<pi in Ei , *F = i(p 2 m S 2 , (3.13) 

where <f>i, <p 2 are constant hermitean matrices (we assume U(N) gauge group here). 
Stokes' theorem gives 

iX=iA=i F = iAifa = - F = -iA 2 (f) 2 . (3.14) 

J C J Si J S2 

The 2d Yang-Mills action localized to constant curvatures is then 

5 = -4" / ^V^tr(*F) 2 = 4- tr (^ 2 + ^(- X ) 2 ) ( 3 - 15 ) 
9 2 d J 9 2 d A i A 2 

and we get the matrix model 

(W) = ^J [dX] ti R e iX e~^ A ^ tlX2 = ^ J [dX}ti R e x e~^ J ^ tlX2 , (3.16) 

where in the last step we have used the relation (12.81) between 2d and 4d coupling 
constants and formally changed variables from X to iX. This is indeed the expected 
form of the familiar one- matrix model (12.10j) . In this derivation one assumes that after 
abelianization the gauge bundle is trivial, which precisely corresponds to dropping the 



contribution of the unstable instantons 43, |49|. In other words, the matrix model 
derived in this way computes the 



perturbative contribution of the 2d YM observables 
as required by the conjecture of {5, @, @]. 

The argument readily generalizes to the case of k non self-interesecting Wilson 
loops of arbitrary topology a general Riemann surface E. We can always dissect E 
into several pieces {E m } which have the topology of a sphere with one, two or three 
holes, such that all given non-intersecting Wilson loops run along some of the cuts. 
The localization to constant curvature gives the multi-matrix model action 

S m . m . = i£-^ tr ( £ s^x) , (3.17) 
^ {Em} ^ m y ieaSm J 

where sf 1 ^ = ±1 is fixed by the orienation of the Wilson loops relative to E m . The 
Wilson loop correlator is then 

(W Rl (d) • • • W Rk (C fc )) = I J [dX x ] ■ ■ ■ [dX k ] tr Rl e iXl • • ■ tr^ e lX * e^ 5 — . (3.18) 
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To compute mixed correlators of Wilson and local operators, we proceed as follows. We 
shrink the Wilson loop Cj which we want to replace by a local operator, and substitute 
in the above action iXi = a,i(*Fi), where is the small area of the corresponding 
disk. We then compute the propagators by inverting the kinetic operator, and take 
the limit of tij — > while killing the propagator from *Fi to *Fi itself to implement 
normal ordering. Inverting again the propagator matrix we then obtain the effective 
matrix model action for mixed correlators with normal ordering prescription in effect on 
local operators. It is easy to verify that this procedure agrees with the matrix models 
shown earlier for the special cases (WO), (WOO) which were derived by summing up 
light-cone perturbation theory. 

While the localization to constant F argument is rather intuitive and convenient, 
it is not completely rigorous. For example, one may worry about subtleties in applying 
Stokes' theorem in the non-abelian case. For this reason, we explicitly show here 
that the matrix model action f 13 . 1 T|) indeed corresponds to the perturbative, or zero- 



instanton sector, contribution to the exact 2d YM partition function 47H51 



To be completely general, let us consider 2d YM with arbitrary compact gauge 
group G on a Riemann surface E (the specialization to G — U(N) is straightforward). 
In the temporal gauge A = the Lagrangian is -^Af, then the Hamiltonian is |g 2 A 
where A is the Laplacian equal to the second Casimir. The exact partition function 
with possibly several Wilson loops inserted is obtained by gluing the spheres with one, 



two or three holes along the boundaries [50J, |51| 



K 1 (U 1 ) = J2dxe-^ AC2{X) Xx(U 1 ) 
x 

K 2 (U U U 2 ) = J2 e ' i92AC2W Xx(Ui)xx(U2) 
x 

K s (Ui,U 2 ,U s ) = J2^e-^ Ac ^xx(U 1 )xx(U 2 ) X x(U 3 ) 
x dx 



(3.19) 



where Ui are the holonomies of the connection on the boundary, A labels irreducible 
representations of G of dimension d\, and xx is the corresponding character. C 2 (X) 
denotes the value of the second casimir of the Lie algebra of G in the irrep A, and A is the 
2d area of a given piece. For example, consider a sphere with three Wilson loops of not 
concentric topology. We can glue this configuration from a sphere with three holes K3 
and three disks K\. The Wilson loop correlator is obtained by including the insertions 
of XR % (Ui) for each Wilson loop and integrating over the boundary holonomies. Explicit 
examples for gauge group U(l) and U(2) are given in appendix. 

We now want to show that the matrix model with action (13.171) is the perturbative 
contribution to the exact partition function obtained from the blocks (I3.19p . For ex- 
ample, consider the block K 3 (the other cases K 2 ,Ki are treated similarly). Let the 
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holonomies Ui be parametrized by elements Xi in the Cartan of the Lie algebra of G, 
such that Ui is a conjugate of an element e Xi of the maximal torus of G. The character 
X\ of the irrep with highest weight A is given by Weyl's formula 

xx(x) = (3.20) 

where W is the Weyl group, R(x) is the Weyl denominator 

R(x) = Y[{e^ a ' x) - e~^ a ' x) ), (3.21) 

a>0 

and p is the Weyl special vector defined by the condition that (p, ctj) = 1 for all simple 
roots0 The dimension of the irrep with highest weight A is 

TT (Ot, A + p) , „ 

*=nH^/- (3 - 22) 

where a > denotes positive roots. The second Casimir C 2 (A) is 

C 2 (X) = (X + p) 2 -p 2 . (3.23) 
Let A(x) be the Weyl measure on the Cartan of the Lie algebra of G 

A(x) = Y[{a,x). (3.24) 

Hence, up to an irrelevant constant factor, the building block K 3 is given by 

Ai^n ^ i d (3.25) 

where P + denotes the subset of dominant weights in the weight lattice of G. 

The exact 2d YM partition function is the integral over the holonomies associated 
with the loops that join the building blocks. Each loop variable enters the measure of 
the integration with the two building blocks associated with it 

Z = J[dU 1 }K Xm (U 1 ,...,)Kx n (U 1 ,...,)... (3.26) 

where S m , S n are two building blocks (either disk, cylinder or triunion) which share a 
boundary loop. We can integrate over adjoint orbits of G and reduce the integral to 
the maximal torus of G on which the measure of integration is 

[dUi] = dxiR(xi) 2 . (3.27) 



11 In finite-dimensional Lie algebras p — \ S Q >o a - 
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The factors R(xi) 2 in the measure cancel the denominators R(xi) in the blocks (I3.25p . 
Let dxi be the flat abelian measure of integration and K &at be the blocks f l3.25|) without 
denominators, i.e. 



(3.28) 

so that the partition function is 

Z = f dxiK™(xi, )K^( Xi , ...,).... (3.29) 

Let us now go back to the matrix model (13.17P obtained from the localization to 
constant F argument. By introducing an auxiliary integration variable $s m for each 
block E m , we can rewrite (I3.17P as 



Here (, ) denotes the positive bilinear form on the Lie algebra g of G which replaces 'tr' 
in the general case. It is normalized such that the long roots a have norm (a, a) = 2. 
We have also omitted the orientation factors sf"^ = ±1 for simplicity!^! While X{ G g, 
it is natural to assume that $s m takes values in the dual g* to the Lie algebra g. To 
proceed, we can integrate over adjoint G-orbits for each variable Xi G g and $s m G g* 
using Itzykson-Zuber formula so that the matrix model integral over X^s and $s m 's 
reduces to the Cartan f) and its dual h*. 

Let X be in the G-orbit of Xi G h, and let $ be in the G-orbit of A G h*. The 
Itzykson-Zuber formula for X G g, $ G g* and integration Dg over the group G with 
the invariant measure reads 

/ Dgeb- 1 ***) = ^ w } t iLi- (3.31) 

JgeG A(x)A(A) 

The matrix model measure j[dXj\ reduces to J dxi with the measure 

[dXi] = J dxiA(xi) 2 (3.32) 

Since each Xi appears in two blocks touching the same boundary, the factors A(xi) 2 
in the measure cancel the denominators A(xi) in the Cartan blocks (I3.3ip . The fac- 
tor A(As m ) 2 in the measure coming from /[d$£ m ] combines with the denominators 



12 We assume canonical orientation of the boundary 9S m with respect to S m , and if necessary 
reabsorb the signs due to orientation into the Wilson loop insertions tr^ ; e ±Xi . 
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A(As m ) in f)3.3ip and gives the overall power A(As m ) 2 ~' l9Sm ', that is exponents 1, 0, — 1 
for disk, cylinder and triunion respectively. Hence, after reduction to the abelian inte- 
grations, the matrix model partition function with action (I3.30P (with possibly Wilson 
loops inserted) can be written in terms of the blocks 

^flat,m. m . = J^dA fl{l2 (-ire< A ^» j j e-^ 2 (3.33) 

and similarly for disk and cylinder. Notice that under the substitution A = A + p, the 
matrix model blocks f !3.33|) are exactly the same as the exact 2d YM blocks (I3.25p . 
except that the summation over the subset of dominant weights P + C ()* in the weight 
lattice P C f)* is replaced by the integration over ()*. Using arguments based on the 
Poisson resummation or Euler-Maclaurin formula one can see that indeed the matrix 
model block (I3.33p . given by the integral, is the perturbative approximation to the 
exact 2d YM block (13 .25 p . given by the sum. To support this general proof, in appendix 
we have computed (PUPFW 7 ) from the matrix model and compared the result to the 
zero-instanton sector of the exact partition function for gauge group U(l) and U(2). 



4 Gaussian three-matrix model in the planar limit 

4.1 Catalan numbers 

For start, we recall the usual combinatorial computation for the generating function 
of Catalan numbers, which count the number of planar diagrams for the one-point 
correlation function tr (X k ^ in the Gaussian matrix model. By c(k) we denote the 
number of such diagrams, so that c(0) = 1, c(l) = 0, c(2) = 1, . . . . The numbers c(k) 
satisfy the recursion on fig. [3] 

k 

c{k) = ^c{i-2)c(k-i). (4.1) 

i=2 

Let f(z) be the generating function 

oo 

f(z) = Y,c(k)z- k -\ (4.2) 

k=0 

From the recursion relation (14.11) . the generating function f(z) satisfies 

f(z) = zf(z)-c (4.3) 
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Figure 3: In this recursion we sum over the label i with which the leg "1" is contracted. 



with Co = 1. Hence the solution with asymptotics f(z) — z + . . . as z — > oo is 



m = i^LH. (44) 



Of course, the generating function 



f(z; A) = \-Wf{\-V*z) = - V f x 4A (4.5) 



is actually the planar resolvent 

N z-X 



/(--A) = (4tr^^) (1.0) 

planar 



in the Gaussian matrix model 

r n 

Z= DXexp( r-trX 2 ). (4.7) 

J 2A 

4.2 One-point resolvent with external legs 

First we compute the number of planar diagrams of the following shape (figHJ). Consider 
a half-stripe / x IR + where an interval I is vertical and R + is a half-line extending to 
the right, and we are given a set P of i + 1 points points placed on / and labelled 
consecutively from "0" to "i" . 

We define a planar diagram with k external legs as a configuration where k points 
from the set P are connected to the positive infinity by horizontal halflines, moreover, 
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C(i'-l) 



C(i'Jc-l) 



Figure 4: stripe diagrams 



these k points must include the point "0" and "i" , and the remaining i + 1 — k points 
are connected pairwise by the internal propagators such that diagram can be drawn on 
the half-stripe without intersections. Let c stI . ipc (i, k) be the number of such diagrams. 

Then 



C s tripc(«, fc) = ^ CstripctV, fc - l)c(i - l' - 1) 



and 



C s tripc(^, 1) — <W 

Let ty(z)fc be generating function for c stripe (i, k) given by 

oo 

= ^2c stT i pc (i,k)z~ l . 

i=0 

The recursion implies = w(z)k-if(z), and hence 



(4.8) 
(4.9) 

(4.10) 



(4.11) 



Now, let c cyc i(n, fc) be the number of planar diagram on a cylinder S 1 x R + with 
the set P consisting of n cyclicly labelled points ... n — 1 positioned on the circle S 1 
and k external cyclicly labelled legs connecting k points from P with the infinity with 
the condition that the legs are cyclicly oriented. In particular, c cyc i(n, n) = n since 
the are n diagrams obtained by cyclic shifts from the diagram in which point "i" in P 
connects by leg "i" to infinity. 

We can deduce c cyc i(n, k) from c str i pe (i, k) by summing over the position of the leg 
"1" and "fc". More explicitly, suppose that the k external legs are connected to points 
on the circle with labels between "0" and "i" (0 < % < n — 1), with the external 
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legs "1" and u k" corresponding to points "0" and "i" respectively. Then the number 
of diagrams is c stT i pe (i, k)c(n — i — 1). Summing over all possible values of "f , and 
including an overall factor of n by cyclic symmetry due to the arbitrary position of 
external leg "1", we get 

n-l 

c cyc i(n,k) =n^2c BtTipe (i,k)c(n-i-l). (4.12) 

i=0 

Consider the generating function for c cyc \(n, k), that is the generating function of 
the one-point diagrams with k > 1 external legs 

oo 

W(z) k = ^2c cy6L (n,k)z- n - 1 (4.13) 

n=0 

From (I4.12p we get 

W(z) k = -d z (w k (z)f(z)) = -d z f k (z). (4.14) 

As a check, from this generating function we can extract for instance c cyc i(l, 1) = 
1, c cyc i(3, 1) = 3,c cyc i(4, 2) = 8,... which can be seen to correctly count the corre- 
sponding number of planar diagrams. 



4.3 Two-point function 

To count two-point planar diagrams we need to glue two effecive vertices W k (z) and sum 
over k. The product W k (zi)W k (z2) overcounts the number of two point diagrams by 
factor of k because on a cylinder we have cyclic symmetry which shifts labelling of the 
k propagators. Therefore, in the sector with k propagators connecting the two-points 
we get the generating function 

W(z h z 2 ) k = ^W(z 1 ) k W(z 2 ) k . (4.15) 

Now consider the two-matrix Gaussian matrix model with matrices X a , a = 1, 2 and 
the propagators 

((X a ))(X b )f) = \ ab 5l5* (4.16) 

After we rescale by the propagators X ab and sum over all k > 1 we get the two-point 
resolvent 

W(zi, z 2 - \ ab ) = d ZuZ2 £ A n )/(z 2 ; X 22 )) k X k 12 (4.17) 



k 
k =i 



and finally 



W(z u z 2 ; X ab ) = -d zltZ2 log(l - AuA/a) = ( tr — tr — ) (4.18) 

Z\ — Ai Z 2 — A2 
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Open 



Figure 5: In these examples, in the closed diagram fc 12 = 4, /c 23 = 2, fc 13 = 3, in the 
open diagram kyi = 3, ki 3 = 3, ^23 = 0. 



where f a = /(s OJ A fl 



57 and used in 



It is easy to check that this agrees with the resolvent derived 
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23[ to obtain exact predictions for Wilson loops and local 
operators. In particular, Laplace transforming the above resolvent on z 2 and setting 
An = (since we want no propagator from the local operator to itself), one finds the 
following result for the correlator of a Wilson loop and a local operator 

(trX/ 1 tr e S2X2 ) conn = J x Ij^y/X^) . (4.19) 

V v ^22/ 



which upon inserting the appropriate values of A12, A 2 2 gives the result obtained in [23 
see eq. (I5.24p below. 



4.4 Three-point function 

When computing the three-point planar connected function we shall distinguish be- 
tween two possible global topologies of the diagrams. We call them closed and open 
type (see fig. Ej). 

4.4.1 Closed type diagrams 

The only difference with the reasoning used to compute the two-point function, is 
that on three-punctured sphere we do not have cyclic symmetry like on a cylinder. 
Therefore, the total number of three-point diagrams is the product of W(zi)k a : 

3 

W(z u z 2 ,z 3 ) = l[W(z a ) ka (4.20) 

a=l 

where k a b, 1 < a < b < 3,a, 6=1. ..3 is the number of propagators by which operator 
a is connected with operator b, and k\ = ^12 + ^23, ^2 = ^23 + ^12; ^3 = ^13 + ^23- Now we 
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will sum over all sectors and also we include dependence on the propagators. Let X a b 
be the propagator between matrices X a and X b . In the sector with fc 12 = 1(^1 + ^2-^3) 
propagators we get a factor A^ 2 2 , etc. Then the resolvent for closed type planar diagrams 
topology is 

W(z 1 ,z 2 ,z 3 ; X ab ) 

fcl2,fc23,fcl3>l ^21,22,23 ^ ^ (-^12/1/2) (-^23/2/3) 3 (-^13/1/3) 

fel2>l,fc 2 3>l,fcl3>l 

(4.21) 

After the summation we get 

w , . , \ _p, Ai 2 A 2 3Ai3(/i/ 2 /3) 2 

W \Z\ } Z 2 , Z 3 , ^o6jfel2,fc23,fcl3>l — ^21,22,23 / 1 , „ „ w 1 , „ „ s/ 1 , £ £ \ 

U — A12/1J2JU — A 2 3/ 2 /3j(l — A13J1/3J 

(4.22) 

with 

fa = f(z a ',\aa). (4.23) 



4.4.2 Open topology 

Now we consider the open topology diagrams (see fig. [5]). Let &4 2 > 1, A; 2 3 > 1, k 23 = 0. 

This case is different from the triangle topology because we have extra diagrams 
with internal propagators connecting legs % and j of operator "1" such that in between 
i and j there could be a bunch of propagators connecting 0\ to 2 . 

The total number of diagrams at operator 0\ with two bunches consisting of k\ 2 
and fci 3 external legs is then 

m— l m—i—2 

Ccycuini, k 12 , k 13 ) = m^2 ^2 ( n i- 1 - i -j) c stripe(i,k 1 2)cstvipe(j^ k n)c(n 1 -2-i-j) 

i=0 j=0 

(4.24) 

In this sum i = i kl2 —i\ and j = ji 3 —ji where i p denotes the leg in vertex "1" connected 
to the propagator p = 1 . . . k\ 2 in the first bunch, and j p denotes the leg in vertex "1" 
connected to the propagator p = 1 . . . k\ 3 in the second bunch. The generating function 

W(*i)*i2,fci3 = Yl c ^2(ni; ki2, k 13 )z^~ x (4.25) 

is 

W(«l)fci2,fci3 = d*iM3l)fcu^2l)fcia(3*/(2l))) 

which simplifies to 



(4.26) 
(4.27) 



21 



Combining W(zi)k 12 ,ki3 w ith the vertices W(z 2 )k 12 an d W(z3)k X3 from operators 2 
and O3 we get the generating function for diagrams with the open connected topology 
in the sector with k 23 = and k i2 > 1, ^13 > 1: 



00 00 



W(zi, z 2 , z 3 ; X)k 23 =o — d 2 zl d Z2 d zz V] V] - — — f^+kn ljfcwjfcia^fe^fcis 

«12 = 1 «13=1 

(4.28) 

We can sum the series using 

00 00 _ 

E E -x^— r*V = — 1^ ^ (4-29) 



l= t ~t * + i - 1 x - y 1 ~ x 

which can be derived from the series expansion of the integral 



1 a+ % y xy 1-y , . 

(1 — tx){l — ty) x — y 1 — x 



So we finally get 



VK(zi, *2, ^3; A af) ) fc2 3 = o = o zl d Z2 d Z3 —- — - log ■ (4.31) 

Al2j2 — A 13/3 1 — M2J1J2 

4.4.3 Complete three-point planar resolvent 

The complete three-point planar resolvent in the Gaussian three-matrix model is given 
by the sum of the the resolvent for the closed topology (14.221) and three possible sectors 
with the open topology (I4.3ip . The final result for the planar connected three-point 
resolvent is then 

W(zx, z 2 , z 3 ; X ab ) = N ( tr l —— tr ^— tr 



Z\ — X % z 2 — X 2 z 3 — X3 



conn 



A12A23A13 (Ziyi's) 2 



-d zl d Z2 d Z3 (1 _ Al2/l/2)(1 _ a 23 / 2 / 3 )(1 - Aia/i/s 

3 2 a a A12A13/2/3 , 1 — A13/1/3 

\A 2 d Z3 —- — — - log - — — — 
A12/2 — A13/3 1 — A12J1J2 



M2J2 — A13/3 

. r, a2 a A12A23/1/3 , 1 — A 23/2/3 

+ d zl d Z2 d Z3 — — — - log — — 

A12J1 — A 23/3 1 — A12/2/1 

, a a o2 A23A13/1/2 , 1 — A13/1/3 

+ 9« d Z2 d zz — — log — — 

A23J2 — A13J1 1 — A23J2J3 



with f a = Za 4 — , a = 1,2,3. As a test, one can verify that if one sets all 
propagators to be equal, A a b = A, then (14.321) reduces as it should to the known 3- 



point resolvent in the Gaussian one matrix model (see e.g. 58]). As a further test, we 
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have also computed (tr X^ 1 tr X^ 2 tr X 3 3 ) from a direct Feynman diagram calculation 
at some fixed low values of ki,k 2 ,k 3 , and verified that the result agrees with the series 
expansion of (I4.32p . 

Note that ( 14.32)) is the resolvent for a general Gaussian 3-matrix model, and hence, 
once the appropriate X ab are plugged in, it encodes the result for all possible correlators 
(WWW), (WWO), (WOO) and (OOO) in our supersymmetric subsector. The 3-point 
function of local operators is a particularly simple case (and of course is well known by 
independent means, since the O's are chiral primaries), so we derive it here as a simple 
test of our general result for the resolvent. Recall that for three local operators, the 
3-matrix model is such that 

A aa = , X ab = — ^ = Xoo for a^b. (4.33) 



Then, (14. 32ft simplifies to 



W( Zl ,z 2 ,z 3 ;X ab ) = d Zl d A *oo(*i*2 + *2*3 + *i*3-2A 00 , 



{^oo — z i z 2){Xoo ~ z 2 z z){Xoo — z l z 3) 

oo oo min(Ji,J 2 ) j 1+ j 2 _ fc , , 

E\T^ A OQ J\J2\J\ + J2 - ^) 

r Jl+l r J 2 +l r Ji+J2-2fc+l K^V 
Ji=U a =l fc=0 Z l Z 2 z 3 

from which can read off the 3-point function (we omit the position dependence on the 
operators, since the correlator does not depend on it) 

J1+J2+J3 

I I s/ \ \ 

(Oj.OjM 




(4.35) 

for I Ji - J 2 1 < J 3 < J\ + -h , «/i + ^2 + ^3 = even 

and zero otherwise. It is easy to check that this is the expected result for the three-point 
function of chiral primaries of the type we consider herein 



4.4.4 Exact < WOO > correlator and strong coupling limit 

We now use the exact planar result to derive the exact prediction for the (WOj-^O^) 
correlator. It will be first convenient to go to an exponential generating function by 
taking the Laplace transform of the 3 point resolvent on zx,z 2 , z$. Let us define for 
convenience 

W(zi, z 2 , z 3 ; X ab ) = d Zl d Z2 d Z3 w(fi, f 2 , / 3 ; A ab ) 

13 The result in (|4.35[) is for the operators normalized as Oj(x) = tr(.T I $ I + i<£ >4 ) J . The 3-point 
function for operators with unit normalized 2-point function is obtained multiplying (|4.35l) by the 
factor (A) 1 



23 



where the function w(f\, f 2 , fa; X a b) can be easily read off from (I4.32p . Then 

N(tr e s ^ tr e s ^ tr e s ^) COQQ = I dz f z ^\ ^ Mz , d Q Q ( ^ 

J {2m) 6 

(4.36) 

where the contours are large circles enclosing the origin. It is convenient to make the 
change of variables z a = X aa f a + l/f a , d Za = Xaa %-i d f a > so that f( z °) = fa, and the 
integral turns into a contour integral over fx, f 2 , 

iV(tr e ^tre^tre^ )conn = _ I ^^ e ^(AuA+^(A aa ii+^)4-3(A»A+^) x 

/ (2-kiY 

x d fl d h d f3 w(fi, f 2 J 3 ; Kb) 

(4.37) 

where the contours are small circles around the origin. Using the contour integral 
representation of the modified Bessel function 

J n=0 v ' 



one finds 

" [/ A L2 \7 A :2:! Yf A,, \\ i + k){i + j)U + k)x 



N{tre SlXl tre S2X ' 2 tre S3X3 ) conn = ^ 

i,j,k=l 



VA11A22/ V VA22A33/ V VA11A33 

x /i +fc (2siVAii)/i+j(2s 2 v / A22)^+fc(2s3\/A3^) 

+ f—^=\ ( \ ik(i + k + 2j- 2)I i+k+2j _ 2 (2 Sl v / A^)/ i (2 S2 v / A^)4(2 S3 v / A^) 

VVA11A22/ VVA11A33/ 

+ f ^i- ) f 7^= ) ^(^ + A; + 2j-2)/ t (2 Sl v^)W2,-2(2s2^)4(2 S3 v / A 
\V A11A22/ VVA22A33/ 

+ f /t~\ — ^ f /^— .) * * (i + fc + 2j - 2)4(2 Sl V ^)/ J (2 S 2V / A^) W2^2(2s 3V / A 
VVA22A33/ VVA11A33/ 



(4.39) 

So far the result is completely general, since we have not specified the propagators X a b- 
For example, plugging in the appropriate values corresponding to the 3-matrix model 
(13.1 Op . the expression f!4.39j) gives the exact prediction for the (WWW) correlator, and 
it is then straightforward to study for example its strong coupling limit. We will not 
pursue this further here, and focus on (WOO) in the following. 

To obtain the exact WO^Oj 2 correlator, we should plug in the A a ;,'s coming from 
( 13~8|) or ( 13~T2|) . and then extract from ffl~39]) (trX/ 1 tr e S2 * 2 trX 3 j2 ) conn , by picking 
up the appropriate power in si,S30 In fact, the result simplifies considerably if we 
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We leave for now an arbitrary S2 parameter. The singly wound Wilson loop corresponds to s 2 = 1. 
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remember the structure of the particular 3-matrix models we want to solve. They are 
such that 

An = A 33 = (4.40) 

since we do not allow propagators from a local operator to itself. Then taking the limit 
An, A33 — » of (14 .391) and extracting the term proportional to s^s^ 2 we get 



N{trX^ ire* 2 * 2 trX^) conn = J X J 2 



min(Ji,J 2 )-l ^Ji-fcxJa-fcxfe 



k=l 
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A 2 



+ 0(J! - J 2 - l)^^|i(Ji - J 2 )/ Jl _ j2 (2 S2 ^) 
A22 2 

00 \J1\J2 

+ E ^rf ( J i + J 2 + 2A; - 2)/ Jl+j2+2fc _ 2 (2 S2V ^ 



fe=i A 



-h+J-2 
22 ' 



+ 6(J 2 - Ji - 1) A V5' ( J 2 - ■/i)/j 1 -j 9 (2s 2 \/A^) 

A 2 2 2 

(4.41) 

where Q(x) is the unit step function (with the convention 0(0) = 1). We can simplify 
this further into the following form 



(trXj^tr e S2X2 trX^ 2 ) 



<nHH yy JlJ2 



A 



12 



22 



A 



23 



min(Ji,J 2 ) 



22 

oo 



X 



fc=l 



where we have defined 



E ^( J i + J 2 - 2k)I Jl+J ^ 2k (2s 2 ^ 2 ) + E(^i + ^2 + 2k - 2)I Jl+ j 2+2k ^ 2 (2s 2 

k=l 

(4.42) 

_ A 22 Ai3 
Ai 2 A 2 3 

Eq. (14.42p is the final result for the correlator {WOO)u computed from the Gaussian 
3-matrix model. To specialize to the two possible topologies of Fig. [2j one should 
insert the appropriate values of X a b given in (I3.7p - fl3.lip . i.e. 



A12 = iX- 



A 23 = ~iX 



A 1 



A 2 ' A 2 

for the open chain OWO topology and 



A13 



A 

A 2 " 



A22 — A 



MM 
A 2 



A\ 2 — — A 23 ; 

for the closed chain OOW topology. 



A 
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A 2 



A22 — A 



MM 
A 2 



(4.43) 



(4.44) 
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We can now study various strong coupling limits of this result. Let us first look at 
the limit 

A -)> oo , J a , J 2 < V\ (4.45) 

Assuming that J\ + J 2 is even (if it is odd, a similar analysis goes through with minor 
modifications), we may rewrite the second sum in (14.42p as (we set s 2 = 1) 

00 00 a 



fc=l 



fc=l 



fc=l 



^(Ji + J 2 + 2A;-2)/ Jl+j2+2fc _ 2 (2v^) = ^(2^(2^)- Yl (2k)I 2k (2^ 2 ) 

(4.46) 



J1+J2 1 
2 

A 22 /i(2 V / A 2 ;)- ^ (2A;)/ 2fc (2 v / A 2 ;) 



fc=i 



Plugging in the contour integral representation (I4.38P of the Bessel function, performing 
the finite sum and evaluating the large A behavior by saddle point analysis, one finds 
that 



£i±£a_! 



A 00 (4.47) 



fc=l - V -*"V^2 2 

In particular, we see that the second term in (I4.46P is subleading compared to the first 
one (recall that I\(x) ~ e x /v2skx at large x). By a similar analysis, one can see that 
the first sum in (14.421) is of the same order as (j4.47p . so the leading behavior of the 
correlator in the limit (I4.45P is 



(trX/Ure^trA^ 2 ) 



J 2 \ 



1 

N 



J\J 2 



M2 \ I A 2 3 



22 



/A 



22 



D 2\A22 



\J 47r\/A 



(4.48) 
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Normalizing by the Wilson loop expectation value 

N 



(tre* 2 ) 



we find 

(trX/Hre^trXf^ 
(tre Xa ) 



1 

iV2 



/A 



/1 2VA 



(4.49) 



22 



A 22 Ji 



A 



12 



Jl 



A 22 J 2 



A 



23 



\ 2 o 



(trX/Ure^)^ (tr e* 2 tr X 3 J2 ) C 



(tre X2 ) 



(tre Xa ) 



(4.50) 



where in the last step we have used the two-point function (I4.19p . We will see in the 
next section that this factorized structure is indeed reproduced in string theory. 
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We can now move on to the more interesting strong coupling limit 

A^oo, Ji/VX = fixed = J x , J 2 < (4.51) 
For comparison to string theory, it is convenient to look at the normalized correlator 

(trX*tre x * trX 3 j2 ) conn 



(trX/Hr e x *) c 



(4.52) 



where the two point function at the denominator is given in (14.191) . To extract the 
strong coupling limit, we can employ as usual the contour integral representation of 
the modified Bessel function and apply the saddle point analysis. The term in square 
bracket in (14.421) gives, after performing the summation and keeping the leading terms 
in the J\ ^> J 2 limit 

^ Up 1 ^ + Ji:^) (4.53) 

2m \ p — z z 1 — z z J 

Writing (see eq. (033]) , (02HD ) 

Jx = Vlji , A 22 = a 2 A/4 , a 2 = 4AiA 2 /A 2 
the saddle point is the solution of 



Ji- + ~(l-^) = 0. (4.54) 

z 2 z z 



The solution which gives the dominant exponential at large A is 



a 2 a 



(4.55) 



Evaluating ( 14.53!) on this saddle point and dividing by the (Oj-AW) correlator we then 
obtain (note that the exponential and the "1-loop" correction around the saddle point 
cancel out in the normalization by the two point function) 

To compare to string theory, it will be convenient to insert also the normalization 
factor ( I2.13P for the local operators, i.e. A/j = ^— i^J=^j ^j- Using the explicit \ a b in 
(I4.43p ~( l4.44[) . we get for the open chain (for which one has p = 1) 



(trX/'tre^trX/ 2 )^. , 1 , . r— ( A, 
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and for the closed chain (for which p = — j^) 



(trXftre x *Af j2 trXf) conn 1 h ( A 2 \ J > /2 



a 1 y a 2 ) 



z~ j2 - Z J J 



+ 



(4.58) 

These are the predictions of the 3- matrix model in the limit (I4.5ip . In the next section 
we will compare them to the result one obtains from semiclassical strings in AdS§ x S 5 . 

Finally, let us look at the limit in which all operators are in the semiclassical regime 

A -»■ oo, Ji/VX = fixed = J x , J 2 /V\ = fixed = J 2 . (4.59) 

Let us assume without loss of generality that J\ > J 2 . Passing to the contour integral 
representation of the Bessel functions and performing the sums in (14 .42 p . one gets the 
following structure 



N(trX^ tr e* 2 trX 3 j2 > conn = j ^- [fr^p^z^^ + f 2 (z)z^ +J ^ 

(4.60) 

where fi(z), fi{z) have at most a linear dependence in VX The saddle point analysis 
at large A then gives the behavior 



1+ (^-J 2 ) 2 _^J2 \ +j2 



iy(trX 1 Jl tre X2 trX 3 j2 ) conn ~ Cl e V V v 

+c 2 e V V V ",(4.61) 

We see that there are two saddle points, both with positive exponents, contributing 
at large A. It can be seen that the saddle point in the first line is the dominant one 
(recall that we are assuming J\ > J 2 ). These saddle points should correspond to 
semiclassical string solutions describing a Wilson loop worldsheet with two "fattened" 
spikes. It would be interesting to find those solutions. 



5 String theory computations 

5.1 Semiclassical < WO > correlator 

For finite R-charge J, it is well known that the correlator (WOj) at strong coupling 



can be computed (as originally done in 59J) via the exchange of a supergravity mode 



between the string worldsheet dual to the Wilson loop and the operator insertion point 
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at the AdS boundary. Here we are interested instead in the limit of large R-charge 
J ~ \/A ^> 1. In this limit the correlator is effectively described by a new string solution 
where the propagator for the supergravity mode has "fattened" into an infinite spike 
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In [26|, the string solution corresponding to (W 1/2 _ BFS Oj) in the large R-charge limit 
was derived (W 1/2 _ BPS is the well-known 1/2-BPS circular Wilson loop). Here we obtain 
a generalization of that solution to the lower supersymmetric case {W 1/4 _ BPS Oj) . The 



Wilson loop Wj/^bps is the quarter-BPS operator of [45j (which is also a special case 
of the infinite family of loops on S 2 of [6]). It corresponds to a circular contour in 
space-time and couples to 3 scalar fields which we take to be $i,$2,^3- The chiral 
primary we consider is given by tr($ 3 + i$ 4 ) J . 

In the original definition of 45|, the 1/4-BPS Wilson loop is specified by the circular 
contour (which we can take to have say unit radius) x^{s) = (cos s, sins, 0, 0), and 
couples to the scalar field combination 

— cos 6*o cos s $i — cos^q sins $ 2 + sin 8 $ 3 . (5.1) 

The choice 9q — tt/2 corresponds to the 1/2-BPS Wilson loop. By a conformal trans- 
formation, one can map the space-time contour x M (s) to be a latitude on S 2 at angle 
6*o : this gives an equivalent definition of Wi/^bps which fits inside the infinite family of 
loops of [6j . To derive the string solution below, it will be more convenient to consider 



the original definition of 45| and do the conformal transformation to the S 2 setup in 
a second step. 

The string solution of interest lies on a AdS^ x S 3 subspace of AdS$ x S 5 , on which 
we take the metric 

ds 2 = — (dZ 2 + dR 2 + R 2 d(j) 2 ) + d$ 2 + sin 2 My 2 + cos 2 Mi\) 2 . (5.2) 

Defining embedding coordinates on S* 3 as 

Q A = (— sin^cosy?, — sin^sinyj, cos $ cos?/', cos d sin ijj), (5.3) 

the local operator we wish to insert corresponds to the spherical harmonic (0 3 + z0 4 ) J = 
(cos$) J e iJ ^. 

We will assume worldsheet conformal gauge, with worldsheet coordinates (r, a), 
and take the following rotationally symmetric ansatz 

Z = Z(r), R = R(t), # = #{t) ^ = V(t) 

(5.4) 

= if = a < cr < 2ir 0<r<oo. 
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The string action specialized to this ansatz reads 



( Z 2 + R 2 + R 2 ) + d 2 + sin 2 + cos 2 ^ 2 



(5.5) 



where the dot indicates r-derivative. 



The solution dual to the 1/4-BPS Wilson loop, without local operator insertion, is 
given by 

Z = tanh r R 
1 



1 



sin $ 



coshr 



(5.6) 



tanh T = sin 9 



cosh(r ± To) 

In these coordinates, the Wilson loop at the boundary sits at r 



0. The choice of 



± sign correspond to the existence of two inequivalent solutions, one stable and one 



unstable, which are dual to the same Wilson loop 45] (see also [14j). For simplicity in 
the following we will focus on the stable solution, corresponding to the choice of +- sign, 
but everything below has a counterpart for the corresponding unstable solution. 

Let us now derive the solution corresponding to the insertion of the local operator 
tr($ 3 + i$4) J , with J/y/X = fixed = J. The insertion of the local operator deforms 
the above solution by creating an infinite spike which terminates at r — > oo. The angle 
ip also assumes a non-trivial profile. The boundary conditions for the various fields are 
as follows 26J. At r = 0, they are controlled by the Wilson loop, i.e. 



7T 



Z(0) = 0, R(Q) = 1, tf(O) = --0 o , V(0) = 0. (5.7) 
At r — > oo, the boundary condition are dictated by the presence of the local operator 



and are given by 26 



]ogZ(r)-tjT, R(t)^0, 0(t)-)-O, ip(r)^ijr, at r^oo. (5.8) 

In particular the boundary condition on ip corresponds to the fact that we have rotation 
on a great circle of S 5 with the appropriate angular momentum J (after analytic 
continuation to Minkowski signature on the worldsheet). 

The equations of motion following from the above action 



-I 







(5.9) 



^ - sin$cos$(l -i) 2 ) = 0, <9 r (cos 2 $V>) = 0, 
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can be solved because they admit the following integrals of motion 



Z 2 1 / • 2 2 
RR Z 

~Z^ + Z~ K (5.10) 

cos 2 flip = p^ 

2 

COS 2 V 

where e, k, p^ are constants. The first and last integrals of motion are just the "ener- 
gies" corresponding to the AdS 3 and S 3 parts of the lagrangian, while the second one 
follows from invariance under dilatation. The fact that the constants in the first and 
last line have equal and opposite values is a consequence of the Virasoro constraints. 
One can see that the boundary conditions in (15.81) further fix p^ = ij, e = J 2 , k = J . 

To integrate the first two lines of (I5.10p . it is convenient to introduce the variables 
(see also |60J) 

R 

u=—, v = log VR 2 + Z 2 (5.11) 

Zj 

in terms of which we have 

ii 2 - u A - (1 + e)u 2 = e - k 2 , v = - . (5.12) 

1 4~ u 

These equations, as well as the last two lines of (I5.1(jp . can be integrated in a straight- 
forward way in terms of elliptic integrals. In fact, the solution becomes elementary 
once one imposes that e = k. 2 = —p^ = J 2 . Going back to the R, Z variables, the 
solution with the correct boundary conditions turns out to be 

Z = r» (5.13) 
cosh(Vl + J 2 (t + T )) 

* %osh(^r + , o) r s ^v^,o)^ 

sintf = , abhCVITjM) = (5.14) 

f\ I Jf J_ 

_ e ~Jr V sin 2 0o sin6» 



1 + -^ 3- 

cos 1? cos i? 



One can see that for 8 = ir/2 (i.e. T\ = oo) one recovers the solution of Zarembo |26 
while for J = we recover the solution dual to the 1/4-BPS Wilson loop (15. 61) . A 
depiction of the AdSs part of the solution is given in Figure El 
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Figure 6: AdS 3 part of the string solution (j5.14j) describing the correlator 
(WOj) for large R-charge J ~ a/A. The infinite spike corresponds to the 
local operator insertion. 



It is not difficult to write down the conformal transformation that maps the solution 
found above to the setup in which the Wilson and local operators are inserted on a S 2 
subspace of the boundary js] 23]. Consider AdS 4 with metric 

ds 2 = -J (dz 2 + dr 2 + r 2 d<p 2 + dxf) . (5.15) 

The transformed string solution lies on the AdS$ slice defined by z 2 + r 2 + x 2 = 1, and 
is given by 

2c Q Z 



1 + c 2 {Z* + E 2 ) 
2c R 

T ~ l + c 2 (Z 2 + R 2 ) 

-i + c 2 (z 2 + r 2 ) e 

X3= l + c 2 (Z 2 + R 2 ) ' c ° = cot 7 



(5.16) 



where R, Z correspond to the solution in f )5.14p . The angle <fi = a as well as the S* 5 
part of the solution are unchanged. One can see that at r = the solution ends on 
a latitude on S 2 of angle 6 Q . Notice that now r = oo is also a point on the AdS 
boundary, since z(r = oo) = 0. Indeed we see that after the conformal transformation 
the spike at r = oo ends on the north pole of the S 2 (i.e. z — 0, — 1, r = 0), which 
is then interpreted as the insertion point of the local operator in the boundary theory, 
see Figure [3 Recall that in the S 2 setup the interesting local operators have the form 
tr (x t ^ l (x) + i$ 4 ) , so an operator inserted at the north pole indeed corresponds to 
tr ($ 3 (a;) + z$ 4 ) , as used in the derivation of the string solution above. 

To obtain the strong coupling prediction for the correlator (Wg Oj), we need to 
evaluate the string action on the solution, supplemented by a boundary contribution 
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AdS 3 




Figure 7: The string solution describing the semiclassical (WOj) corre- 
lator after the conformal transformation (15.161) . The Wilson loop is a 
latitude on S 2 = dAdS^, here depicted as a plane, and the spike ends on 
the north pole of S 2 , where the local operator is inserted. 



at infinity which follows from the insertion of the local operator at r = oo. After using 
in ( 15. 5p the values of the integrals of motion, the "bulk" action is 



'bulk 



V\Jdr(^ + sin 2 =V\J dr(^ + sin 2 A . (5.17) 



This has a well- understood divergence at r = 0, which in this case we can simply 
eliminate by considering the normalized correlator (Wg Q Oj)/ (Wg ), i.e. by subtracting 
the value of the action at J = 0. After performing the r-integral 



(5.1* 



POO 

/ dr (sin 2 - sin 2 &\ J=0 ) = sin # - 1 + y/l + J 2 - \/ J 2 + sin 2 6 
Jo 

we obtain 

Sbuik - S bu ik\j=o = \/A (sin 6»o - V J 2 + sin 2 6U . (5.19) 



The boundary contribution at infinity originates from the insertion of the vertex oper- 
ator for the local operator at the north pole, which, omitting subleading factors in the 
large J = \f\J limit, is essentially given by 



z 2 + r 2 + (x 3 - 1)V v ' \2 
This corresponds to adding to the action the boundary term 



(cos^) J e iJ ^ = (— (cos??) J e^. (5.20) 



Sbdy = -V\(J\ogZ(oo)+iJil;(oa) + J\ogcos#(oo)) - V\J\ogj. (5.21) 
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The last term in the first parenthesis vanishes, while the divergences in the first two 
terms cancel each other leaving the result 

S bdy = -V\J\og(Jl + ) - VXJ log ^ . (5.22) 

\ y sin 6*0 sin I 2 

Putting everything together, we then find 



(W do Oj) = e-^-^y = e v^^^ 2 Weo-sineo+Jiog^i+-^-^jj+VAiog(ico t ^) 
(Wg ) e" 56 «' fc lj=° 

(5.23) 

We can now compare this to the field theory prediction. The relation to 2d YM implies 
that the exact Wilson-local correlator should be given by the Gaussian two-matrix 



model result in eq. (I5.24p . see |23J. Inserting the values of the couplings Ai 2 , A 22 which 
can be read off from the action in ( 12.111) . as well as the normalization factor (12.131) for 
the local operator, (I4.19P yields the prediction, for an arbitrary Wilson loop on S 2 and 
local operator inserted at the north pole (or any other point on the upper hemisphere) 



(Oj(x N )W(C)} _ 1 ^ / ^A * Ij(VX') (5 24) 



(W(C)) 2 \AAJ /.(VA 7 ) 

where A' = \A X A 2 IA 2 \. Recall that for a latitude AAiA 2 /A 2 = sin 2 9 and A 2 /A 1 = 
cot 2 f . Thus we see that the prefactor {^j = e ^J^\^MlM = e vXJio g (±cot^) 
reproduces the last term in the exponent in (I5.23p . The remaining terms arise from 
the strong coupling limit of the Bessel function 

sin 9 ) = I *L z J-i e ^**M*+±) . (5.25) 

The saddle point analysis at large A and large R-charge J = J[y \ is the same as in 
(I4.54p - p.55p . and the dominant exponent at strong coupling^ gives 

Jj(VAsin0 o ) ~ e v ^( v/j2+sin2eo+Jlog ( v/l+j2/sin2eo " J/sin90 )) . (5.26) 

After dividing by ii(-\/Asin#o) ~ e ^ sm6 o^ we ^ ere f ore precisely reproduce the string 
theory prediction ( 15.231) . 



5.2 < WOO > correlator 

In this section, we present the string theory computation of the correlator (Wg Oj 1 Oj 2 ), 
and we then compare to strong coupling limit of the 3-matrix model prediction. 



15 There is also a subleading saddle point which corresponds to an unstable string solution, see refs 
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Figure 8: The two diagrams contributing at leading order at large N to 
(W O j x O j 2 ) for Jx, J 2 <C vA. Fig. (b) is subleading in the large A limit. 



Let us start with the case in which both local operators are "light" , i.e. J\, J 2 <C y/X. 
In this regime, the correlator at strong coupling should be computed by considering 
the exchange of two supergravity modes between the boundary and the worldsheet 
describing the 1/4-BPS Wilson loop 15.61 There are two types of contributions, one 
involving a 3-point interaction in the bulk, and one with two independent propagators 



ending on distinct points on the worldsheet, see Figure El As explained in |44j in the 
case of the analogous calculation of 4-point functions of local operators, the diagram 
involving the 3- vertex is subleading at large A, since it involves only one vertex operator 
inserted on the worldsheet. On the other hand, the diagram with two independent 
bulk-to-boundary propagators simply factorizes, so we immediately conclude that 

(5.27) 

Each of the factors on the right-hand-side is the same as computed in [23[, in particular 
it agrees with the two-point function in the two-matrix model. The factorized structure 
precisely agrees with the 3-matrix model in the relevant limit, see eq. (I4.50P 1^1 It would 
be interesting to compute the subleading corrections from the diagram in Fig. 8(b), 
and compare to the corrections to (14.50P on the matrix model side] 17 ! 

Next we look at the more interesting case in which one operator is "heavy" J\ ~ 
vA 3> 1, while the other is "light", i.e. J 2 is kept of order one. In this case, we can 
simply use the string solution for (Wg Oj x ) derived earlier, and compute the amplitude 
for the exchange of a supergravity mode corresponding to Oj 2 , see Figure EB The recipe 



to do this is well understood, see e.g. 30[ 32( for analogous calculations involving 



three-point functions of local operators. Here we follow the notations of [23J, where 



16 It is easy to see that a similar leading factorization into a product of (WO) applies more genereally 
to the correlator of a Wilson loop and several local operators, as long as the operators are "light" . 



17 



For a recent evaluation of the analogous subleading corrections to the 4-point function of two 



•'heavy" and two "light" operators, see [61 1 
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Figure 9: Diagram describing the string theory computation of the cor- 
relator (WO^O j 2 ) for Ji ~ y/\ and J 2 <C J\. 



the appropriate vertex operator for local operators inserted on S 2 was given. We have 
(W 6o Ojtf N )O j2 {2o)) 



j d 2 a (y AdS ^ + V 



.1-2 



Vg* = ^ ^t- 1 ( V 2( / 2 y ^d^d^Gj^Xo) - -^d a d a G, h {x»-x Q ) ) Y J %e,x ) 



T/5 5 Vx j 2 + i 
Q A 



■ 2 V J 2 + l ^ 2 ' J 2 + l 

2 J 2 d a Q A d a Q A G j2 (x»; x )Y J *(O, x ) 



G 



A-nN AyfJ 2 
c) , £ = ( 

sin 1? cos <y9, — sin 1? sin ip, cos cos ifj, cos 1? sin ^) 

Jn / _ \ -/a 

yJ 2 



2 2 + — Xo) 



2 — 2x ■ x~o 



x -B + i& 



■h 



(5.28) 

where A correspond to the solution (15.141) . with z, r, X3 related to R, Z in ( 15. 14ft by 
the conformal transformation (I5.16P (and of course J = J\ = Ji/\/\). The notation 
xn indicates the north pole of S 2 , where the heavy operator is inserted, while we 
keep the insertion point of the light operator xq G S 2 arbitrary. We can choose the 
parameterization 



X 



(sin 77 cos 7, sin 7] sin 7, cos rj) 



0<r/<7r, < 7 < 2vr. 



By rotational invariance of the problem, we can always set 7 = 0, which we will assume 
in the following. We can simplify these expressions further if we integrate by parts the 
second term in V^ dSs . This produces a term proportional to d 2 Y J ' 2 . Using the classical 
equations of motion d 2 Q A + (d a Q B d a Q B )Q A = 0, we can rewrite this as 



d n d n Y J z 



J2O/2-1) (d a hg(x -Q + ie 4 )) -j 2 d a e A d a e 



Y J 2 



(5.29) 



Then, further using the fact that on the explicit solution -^daX^daX^ = 2^ + J 2 and 
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d a Q A d a Q A = 2 sin 2 & — J?, we arrive at 
(W eo O. h (x N )O j2 (x )) 



(Wg Oj^Xn)) 



4tiN 
'^NJ 2 J 



(J2 — 1) / drda 



r l f / -> \ \ 2 



— -sin^ + J"/- o log fo-0 + *© 



^2 



T=0 



(5.30) 



The boundary term in the second line is easily evaluated. The contribution from r = 0, 
where the Wilson loop is inserted, vanishescJ On the other hand at r — > 00 we have 



-JlT 



J2 



G 



■h 



c (l-cosr ? )( v ^T^- Ji) 



+ . . . 



Y 



-(1 — cos7])e JlT 



(5.31) 



'1 + J(- Ji 



7 2 

1 + 



+ ... 



sin 2 0o sin #0 



so that 



2vr 

rfa [(9 T G j2 )F j2 -G j2 (9 T r j2 )] 



r=0 



tanf 



'1 + 



J? 





sin 60 . 

(5.32) 

Note in particular that the dependence on the insertion point Xq of the light operator 
(which enters through the parameter rj, see above) drops out. 

The integral in the first line of (I5.30p is more challenging to compute, given the 
complicated form of the solution. A limit in which it can be evaluated analytically 
is 9o = I and J\ — > (in this case the calculation essentially reduces to the (WO) 
correlator considered in 23]). In particular, in this limit one can explicitly show the 
independence on the position of the local operator Oj 2 , see appendix B. In any case, 
a numerical evaluation is straightforward for all ranges of the parameters. Keeping 
in mind the structure of the matrix model result fl4.57l) - fl4.58l) which we would like to 
match, we can simply guess the analytic result and compare it to numerics. We find 
that the numerical integration agrees with very good precision, over a wide range of 
the parameters J2, Ji, , rj, with the following analytic expression 



4nN 



(J2 



drdo 



- sin 2 ?? + Jl - (<9 a log (x - © + i© 



Gj 2 Y J * 



18 Recall that J? > 2. 
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2 72 v^A(-tan|) J2 Ji^J ^o<^<tt 

(cot *)* * p--^^ + a) »<„< r> 

Here we have used, as in the matrix model calculation, the shorthand notation 



Ji Ji 



sin 9 srn #o 

Note that, as a non-trivial result of the integration, we find that the correlator is 
independent of the insertion point Xq, except for a jump when the operator crosses the 
loop. This is again as expected from the localization to the 2d theory. 

Putting together the integral (I5.33P and the contribution of the boundary term 
(I5.32p . our final result for the 3-point correlator from string theory is 

(W 9o Oj 1 {x N )O j2 {x Q )) 



-J 2 



^v^Hanf/ 2 Jif£ x eS- 

N 



iv^A (cot f ) J2 J, tan 2 f X^*_ - + ^ * G \ fe} 



where 5 denote the two regions of 5" 2 singled out by the loop, with the north pole 
xn G S + . Recall that for a loop along a latitude with angle 9 we have 



A 1 = 4 7 rsin 2 y, A 2 = 4vrcos 2 y, (5.35) 

and so we see that the result of the string calculation precisely matches the 3-matrix 
model prediction fl4.57p - fl4.58p . 



5.3 Supersymmetry equations 

In this subsection we show that the string solution in f)5.14p dual to the (WO) correlator 
preserves the same supercharges found in the field theory analysis. Moreover we derive 
a set of first order differential equations which should be satisfied by the string solutions 
dual to correlators of any number of Wilson loops and local operators. 

19 It might be possible to prove the position independence of the correlator for a generic solution in 
the supersymmetric subsector on S 2 . One should be able to represent the variation of the integrand 
under the shift of the position xq € S 2 of the vertex operator as an exact two- form using the calibration 
conditions (see the next subsection) and the fact that the operator SO, where 5 is a shift along S* 2 , is 
Q-exact. 
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Let us start with the K-symmetry equation for a fundamental string. Since the 
relevant solutions lie on AdS% x S 3 , we can restrict the analysis to such subspace of 
AdS§ x S 5 . It will be actually convenient to view AdS% x S 3 as embedded in AdS^ x S 3 
with metric 

ds 2 = \dx i dx i + z 2 (dy i dy i + (dy 4 ) 2 ) % = 1, 2, 3 
z- 2 ee yY + (y 4 ) 2 . 
The equation defining the relevant AdSs x S 3 subspace is 

X V + z 2 = 1 . (5.37) 

In what follows, we will denote the coordinates on AdS± x S 3 as X M = (x\y\y A ), 
M = l,...,7. 

In our conventions and in Euclidean signature, the ^-symmetry equations read 

id T X M d a X N T MN e AdS = Vhe AdS . (5.38) 

Here T MN = 1/2\T M , T N ] and T M are curved 7d Dirac matrices, e A ds is the AdS 5 x S 5 
Killing spinor and h is the determinant of the worldsheet metric. We assume worldsheet 
conformal gauge, so that 

Vh = d T X M d T X N G MN = d a X M d a X N G mn , (5.39) 

where Gmn is the metric (I5.36p . 

Restricted to AdS^ x S 3 with metric f)5.36p . the explicit form of the Killing spinor 

is 

e AdS = -j= (e + z[x i Y i - y% +3 - y 4 r 7 ) ei ) = (e + X M T M ei) , (5.40) 

where for later convenience we have introduced the notation 

X M = z{x\-y\-y'). (5.41) 

Note that on the AdS% x S 3 subspace we have X m Xm = 1- Note also that X M is 
orthogonal to all vectors in the tangent space of AdS^ x S 3 

X M p M = Vp M e T{AdS 3 x S 3 ) (5.42) 

as follows from differentiating the constraint x l x % + z 2 = 1. 

In ( 15.401) eo and ei are constant 16 components (Majorana-Weyl) spinors which are 
directly identified with the Poincare and superconformal supersymmetry parameters 
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in the M = 4 SYM at the boundary. In particular, to impose that the string solution 
preserves the same supercharges as the dual operators, e and e\ should satisfy 

iliP^o = Piti , iliPiti = Pito 

which imply that 

(Tu + Pij) e o = , 7i2 3 eo = P46 ^ 
ei = -ip^o • 

We can collect all these conditions in the following equations, written in terms of curved 
7d Dirac matrices 

i^MN^o — JmnpF P €i , iTmn^i = Jmnp^ P ^o ■ (5.45) 

Here Jmnp is a totally antisymmetric tensor whose components can be read off from 
fl5.43p - fl5.44p . They are the components of the following 3- form 



J (3) = ±. j MNp dX M A dX N A dX p 



1 

3! 

1 z 3 ■ 

—e ijk dx l A dx j A dy k - —t ijk dy % A dy j A dy k - zdx 1 A dy % A dy A . (5.46) 



Similarly to [6J, from J( 3 ) we can construct an almost-complex structure on AdS^ x S . 
Indeed, one can see that 

J A/ = P X P (5.47) 

satisfies on AdS^ x S* 3 

J A P = -6% + X M X N (5.48) 

and hence it squares to —1 when acting on a vector in the tangent space T(AdS% x S 3 ). 
It is also easy to see that maps tangent vectors to tangent vectors, and so it defines 



an almost-complex structure on AdS% x 5* 3 o From J A ^ one can construct as usual a 
2-form J(2) = l/2JMNdX M A dX N , which however is not closed. 

Going back to the ^-symmetry equations, we start from 

(id T X M d c X N T MN - Vh) (e + X p r pei ) = (5.49) 
and use the supersymmetry conditions ( 15. 45 p . This gives 

(d T x M d a x N x Q v Q j MNP v p - Vh) eo + (uV ^pr p - v^x P r p ) £l = , 

20 In fact, its components are closely related to those of the almost-complex structure on AdS± x S 2 
constructed in 0] for general supersymmetric Wilson loops on S 3 . 
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where we have used that \Tmn,Tq] = IGnqY m — 2GmqTn and that Xm9 t X m = 
X M d a X M = 0. Writing r Q r p = 5q + T Q P and using again (OS]) as well as the last 
line of f)5.44p we end up with 

d T X M d a X N J MN -Vh)l (5.50) 



eo = 0. 



■-3 [drX M d a X N J m n p — VhXp — id T X M d a X N JmnqJ%^ r p r 456 



Since we have already imposed all supersymmetry constraints, the terms multiplying 
the identity and r p r 456 should vanish independently The condition coming from the 
identity 

d T X M d a X N J MN -Vh = (5.51) 

is the statement that the string solution is "calibrated" by J( 2 ) (in a generalized sense, 
since J@) is not closed), i.e. the worldsheet area is given by 

A(E) = f J {2) . (5.52) 



Note that the calibration condition is actually not independent from the remaining 
seven conditions 

d T X M d a X N J M NP - VhXp - id T X M d a X N J MNQ J Q p = (5.53) 

since (I5.5ip follows from contracting the above equations by X p (recall that X m Xm = 
1 and J Q P X P = 0). 

The fact that the string solutions are calibrated by J(2) is analogous to the result 
obtained in for the Wilson loops on S 3 . There is however an important difference 
in the present case involving mixed correlators of Wilson loops and local operators. If 
the string solution X m (t, a) were real, then from 

{J M N d a X N - d T X M f = 2Vh- 2J MN d T X M d a X N = (5.54) 

one could conclude that the solution is pseudo-holomorphic with respect to J, i.e. 

J M N d a X N = d T X M . (5.55) 

One can also check that if (I5.55P holds, then f!5.53j) are automatically satisfied!^! 
So for real solutions, the supersymmetry equations are equivalent to the pseudo- 
holomorphicity equations. This is the case for correlators of Wilson loops alone, as 



21 To see this, one can note that J( 3 ) satisfies the identity JmnqJ % = GmpXn — GnpX m + 

^e MNPQlQ . 2Q3Qi xQijQ^Q*. 
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in |6j. However, in the cases involving some local operators, the boundary conditions 
on the Euclidean worldsheet are such that some angle along S 5 is purely imaginary, see 
([52]) Ej, so the solution X M (r, a) is not in general real. Therefore, while the solutions 
are still calibrated (15.521) . the pseudo-holomorphic equations do not follow from (15.541) . 
and one should solve the more general supersymmetry equations (15.531) . 

To support the general analysis of this subsection, we have explicitly checked by 
direct computation that the solution (15. 14}) found earlier indeed satisfies (15.531) . proving 
that it preserves the correct supersymmetries. On the other hand, it is easy to see that 
it does not satisfy the pseudo-holomorphicity equations (I5.55p . unless J = 0. 

It would be interesting to study further the supersymmetry equations (I5.53P and 
see if they can be used to find new solutions dual to more general (WW ■ ■ ■ 00 • • • ) 
correlators. In particular they could be useful to find the solution dual to (WOj-^O j 2 ) 
for Ji, J 2 ~ V% for which we have an explicit prediction (14.611) from the 3-matrix 
model. We leave this for future work. 
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A < WWW > for U(l) and U(2): matrix model vs 
exact 2d YM partition function 

A.l U(l) 

Let G = U(l). Let the holonomies be Ui = e lWiXi . Let the Wilson loops be e lWiXi where 
Wi are fixed integer weights. 

22 While the analysis in the previous section was for (WO), a similar boundary condition would 
apply also in the general case, at each spike corresponding to a local operator. 
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Gluing K tT i and three i^disk and integrating over X{ we get 

(W 1 W 2 W 3 ) U{1) = — ^ e -\{n 2 a+(n+w 1 fa 1 + {n+w 2 ) 2 a 2 +{n+w A fg. i ) ^ 

where here n labels the summation index over irreps in the block K tTi , and for shortness 
we absorbed g 2 into areas a, where a is the area of the triunion and a, are areas of 
the disks. Now we will make use of the Poisson resummation identity 

J25(y-n) = J2^ mmy , (A.2) 

so that 

<WiW 2 W 3 ) m = ± I dyj^e 271 

1 v-^ (^a W + 27rtm)^ t ^ , 
= —2_^ e a+2>i i^i W i a 



rimy ^-\(y 2 a+(y+w 1 ) 2 ai+(y+w 2 ) 2 a2+{y+w 3 ) 2 a- i ) 

(A.3) 



where now the meaning of m is labelling 2d magnetic fluxes, or 2d unstable instantons. 
The perturbative, or zero-instanton, contribution comes from the m = term: 



<WW 2 WM rU/m = e^^ E ^ (A.4) 



'pert,U(l) 



If we repeat the localization argument of 14J to the constant curvature configurations 
as described in section I3.2[ we get the Gaussian 3-matrix model 

(W 1 W 2 W 3 ) pert = J dX 1 dX 2 dX 3 e~^ ± ^ 1 ~ j: ^W(X l )W(X 2 )W(X 3 ) (A.5) 

where Xi take value in the Lie algebra of G, and Wilson loops are defined on the group 
elements e lXi G G. Not surprisingly, the Gaussian integration of the U(l) matrix model 
with W(xi) = e lWiXi , using the standard formula J e -k xHx + lJx = e - 2 JH lj ; reproduces 
the perturbative part of the correlation function from the exact Migdal formula for 
U(l) (TOJ. 



A.2 SU(2) 

We will label irreps of spin j of SU (2) by their dimension n — 2j + 1, so that n e Z>o- 
The characters are 

, , sin nx , . . 

Xn{x) = A.6 

SIM 



The Casimir is 



C 2 = 2j(j + l)= 1 -(n 2 -l) (A.7) 
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(our Casimir is twice of the standard quantum mechanics textbook conventions £V L\ = 
+ because of the normalization tr F L 1 L 1 = tr F L 2 L 2 = tr F L 3 L 3 — \)- The build- 
ing blocks are 



_n2-i 1 sin nxi sin nx 2 sin nx 3 
K tri (x 1 ,X 2 ,X 3 ) = ) e * : : : 



n smxi smx 2 smx 3 

n=l 

oo , 

jr / \ \- _!!izi a . sinnjXj 
^— ' sm Xj 

n;=l 



(A.8) 



The Haar measure to integrate over G = SU (2) is 

Let the Wilson loops be in the fundamental W = \tv F U(x) ) i.e. 

Wi(xi) = cosxj. (A. 10) 

Integration over Xi is trivial 

J sinnx sinnxicosxi = ^(<5„, ni +i + 5 n , ni -i) ( A -H) 

so that we obtain 

(WiW a Ws> w(2) = -- ^ )^exp ( --a - ^ a, ) -L^- 



'St/(2) - Z g P I g 

(iui,w 2 ,w 3 )=(±l,±l,±l) n=1 V i=1 

= E E e ^ a " n/3_7 ( n2 + nc iH + c 2H + - c 3H) ( A - 12 

o 71 

(uii,?i)2,'i«3)eS n^O 



where 



and 



a = ^(a + ai + a 2 + a 3 ) /3 = ^^u>ia; 7 = ^E w2ai (A.13) 

ci(iy) = c 2 (u>) = WjWj c 3 (w) = Wiw 2 w 3 ( A - 14 ) 

and the index set S = {(1, 1, 1), (—1, 1, 1), (1, —1, 1), (1, 1, —1)}. Next we make a Pois- 
son resummation (there is a subtlety as we have sum over n ^ not over n G Z but 
the extra term n = actually gives vanishing contribution after summation over S). 
We get 

(w 1 w 2 w 3 ) su{2) = \\ - c iH^ + C *H - V) e ^ +WBN 

(u>i,u> 2 ,u> 3 )es ™s2 

(A.15) 
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(W 1 W 2 W 3 ) perttSm = 1 - V — A , A , ; , A / + • • • (A. 16) 



where 9 (3 1 F(/3) = d(3F((3). The perturbative part is given by the term m = 0. We 
find 

3 (A(E Ai) + Ek,- 2^^) 
8 (A + Ai + A 2 + A 3 ) 

If we wish to get the U(2) correlation function, we simply need to multiply the above 
SU(2) result by the U(l) correlation function (1A.4I) with g 2 substituted by \g 2 (be- 
cause in our conventions the Casimir of U(l) factor inside U(N) is ■h in fundamental 
representation). 

Starting from the matrix model ( IA.5j) for U (2) (i.e. the matrix integral is over 2x2 
hermitean matrices), we have computed perturbatively (WfWfWf) up to order g 10 
and verified that it agrees with the expansion of the perturbative truncation HA.16[) 
(multiplied by the £7(1) contribution (IA.4|) as described above) of the exact result 
( 1A.15j) . Of course, this agrees with the general proof given in section [3721 



B Position independence of the < WO > correlator 
from string theory 

Here we show analytically the independence of the integral in f !5.33|) on the position xo 
of the local operator in the special case 9q = | and J\ — > 0. In this case the calculation 



reduces to (WO) and the non-trivial integral to compute is (see also 231 ]) 

1 p2tt roc j / 

T_ 1 f ^ f dT ^U-2-( COST l 



, da 2-tanh 3 -V : =- (B.l) 

2tt Jo Jo cosh r \ 1 — sin rj cos a cosh t J j — 1 

for any integer j > 2. The integral is elementary for rj = by the change of variables 
£ = tanhr. We want to show that, in fact, the integral does not depend on rj. In terms 
of the variables sin {3 = cosh -1 r, cos (3 = tanh r the integral is 

™ = (* *, P v*±L( c03 " co ^. y . (b.2) 

J J Q cos z p \ 1 — cos a sin rj sin p J 
We now notice that the integrand 

/ cos.cos^ y 

cos^ p \ 1 — cos a sin rj sin p y 

satisfies 

cosrj drjf da A d/3 = dg (B.4) 
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where the one-form g is 



g = f cos cr cos (3 da + / -dj5 . (B.5) 

sin p 

Therefore 

cosr]d v (27rI) = dg = g = (B.6) 
Jr J or 

where the cylinder R = Ip x is the integration domain and dR is the boundary 
consisting of two circles Si at (3 = and /3 = | , and the last equality holds because g 
identically vanishes for = or /3 = |. Hence we have shown that J does not depend 
on in the range < rj < | in which /, g are smooth on the compact integration 
domain i?. 
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